Logarithmic and Riesz Equilibrium for Multiple Sources on the Sphere --- the Exceptional Case

We consider the minimal discrete and continuous energy problems on the unit sphere $\mathbb{S}^d$ in the Euclidean space $\mathbb{R}^{d+1}$ in the presence of an external field due to finitely many localized charge distributions on $\mathbb{S}^d$, where the energy arises from the Riesz potential $1/r^s$ ($r$ is the Euclidean distance) for the critical Riesz parameter $s = d - 2$ if $d \geq 3$ and the logarithmic potential $\log(1/r)$ if $d = 2$. Individually, a localized charge distribution is either a point charge or assumed to be rotationally symmetric. The extremal measure solving the continuous external field problem for weak fields is shown to be the uniform measure on the sphere but restricted to the exterior of spherical caps surrounding the localized charge distributions. The radii are determined by the relative strengths of the generating charges. Furthermore, we show that the minimal energy points solving the related discrete external field problem are confined to this support. For $d-2\leq s<d$, we show that for point sources on the sphere, the equilibrium measure has support in the complement of the union of specified spherical caps about the sources. Numerical examples are provided to illustrate our results.Comment: 23 pages, 4 figure

Solving parabolic equations on the unit sphere via Laplace transforms and radial basis functions

We propose a method to construct numerical solutions of parabolic equations on the unit sphere. The time discretization uses Laplace transforms and quadrature. The spatial approximation of the solution employs radial basis functions restricted to the sphere. The method allows us to construct high accuracy numerical solutions in parallel. We establish $L_2$ error estimates for smooth and nonsmooth initial data, and describe some numerical experiments.Comment: 26 pages, 1 figur

Application of quasi-Monte Carlo methods to PDEs with random coefficients -- an overview and tutorial

This article provides a high-level overview of some recent works on the application of quasi-Monte Carlo (QMC) methods to PDEs with random coefficients. It is based on an in-depth survey of a similar title by the same authors, with an accompanying software package which is also briefly discussed here. Embedded in this article is a step-by-step tutorial of the required analysis for the setting known as the uniform case with first order QMC rules. The aim of this article is to provide an easy entry point for QMC experts wanting to start research in this direction and for PDE analysts and practitioners wanting to tap into contemporary QMC theory and methods.Comment: arXiv admin note: text overlap with arXiv:1606.0661

Matching Schur complement approximations for certain saddle-point systems

The solution of many practical problems described by mathematical models requires approximation methods that give rise to linear(ized) systems of equations, solving which will determine the desired approximation. This short contribution describes a particularly effective solution approach for a certain class of so-called saddle-point linear systems which arises in different contexts

Stability and preconditioning for a hybrid approximation on the sphere

This paper proposes a new preconditioning scheme for a linear system with a saddle-point structure arising from a hybrid approximation scheme on the sphere, an approximation scheme that combines (local) spherical radial basis functions and (global) spherical polynomials. In principle the resulting linear system can be preconditioned by the block-diagonal preconditioner of Murphy, Golub and Wathen. Making use of a recently derived inf-sup condition and the Brezzi stability and convergence theorem for this approximation scheme, we show that in this context the Schur complement in the above preconditioner is spectrally equivalent to a certain non-constant diagonal matrix. Numerical experiments with a non-uniform distribution of data points support the theoretically proved quality of the new preconditioner. © 2011 Springer-Verlag